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Abstract 

We present transport calculations for heavy ion reactions in which the mean field 
and the in-medium nucleon-nucleon cross section are consistently based on the same 
effective interaction, i.e. the in-medium T-matrix from microscopic Dirac-Brueckner 
calculations. Doing so, the stopping in central reactions in terms of the recently 
proposed var t \ observable and the correlation to the behavior of the directed flow 
^-j- | is investigated. The relation to the nuclear shear viscosity is discussed. 
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X 1 One of the primary goals in studying heavy ion reactions at intermediate ener- 

gies is the investigation of the nuclear equation-of-state (EOS) at supra-normal 
densities and/or high temperatures [1]. In a hydrodynamical picture the time 
evolution of such a reaction can be understood in terms of pressure gradients 
which build up in the compressed zone and drive the dynamics. Originally it 
was therefore expected that a direct determination of the nuclear EOS via col- 
lective flow observables should be possible [2]. However, the situation turned 
out to be more complex: The system does not behaves like an ideal fluid but 
binary collisions of the nucleons lead to a viscous behavior. Moreover, over 
most of the reaction time, in particular during the compressional phase, the 
system is even out of local equilibrium [3,4]. E.g. experimental evidence for in- 
complete stopping even in central reactions has been reported in [5]. Thus the 
hydrodynamical limit is not reached in relativistic heavy ion reactions except 
of the extreme ultra-relativistic regime [6]. 

On the other hand, microscopic transport models turned out to be an ade- 
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quate tool for the description of the reactions dynamics at intermediate ener- 
gies. The physical input of such Boltzmann type semi-classical models are the 
nuclear mean field U and the nucleon-nucleon (NN) cross section a. Both are 
determined by the effective two-body interaction in the medium, i.e. the in- 
medium T-matrix; U ~ 9?Tp, a ~ QT, respectively da/dVl ~ |T| 2 . However, 
in most practical applications phenomenological mean fields and cross sec- 
tions are used. Adjusting the known bulk properties of nuclear matter around 
the saturation point one tries to constrain the models for supra-normal den- 
sities with the help of heavy ion reactions [7,8]. Medium modifications of the 
NN cross section are often disregarded which works, from a practical point of 
view, astonishingly well. However, in particular kinematics regimes a sensitiv- 
ity of dynamical observables such as in-plane (i>i) and elliptic (1*2) flow [9] or 
transverse energy transfer [10] to the NN cross section has been observed. 

In this letter we present calculations were both, the mean field and the in- 
medium cross sections are consistently based on the relativistic Dirac-Brueckner 
(DB) T-matrix. The standard Dirac-Brueckner ladder approximation [11] ac- 
counts for two body-correlations to lowest order in the Brueckner hole-line 
expansion. Microscopic DB mean fields have already been successfully applied 
to heavy ion reactions [4] but a consistent treatment of mean field and cross 
section has not been performed so far. In a consistent truncation scheme two- 
body correlations are maintained at the same level in the drift term and the 
collision integral of the kinetic equation [12]. 

At finite densities the NN cross section is suppressed for two reasons: (1) 
Kinematics, i.e. the cross section scales with the in-medium nucleon mass as 
(7 e fr ~ (M* /M) 2 a [13,14]. (2) At low energies an additional suppression due to 
Pauli-blo eking of the intermediate states in the Bethe-Salpeter (BS) equation 



appears. The Pauli-Operator Q — (1 — /)(1 — /) in (1) projects intermediate 
states onto unoccupied phase space regions. The mean field, respectively the 
Hartree-Fock self-energy, is obtained by a summation over the Fermi-sea S = 
— itr[TGp]- This leads to dressed quasi-particles carrying effective masses and 
kinetic momenta 



given by the real part of the, in general complex, self-energy. For on-shell 
scattering the NN cross section is given by 



T — V + iVGGQT 
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Fig. 1. Elastic in-medium neutron-proton cross sections at various Fermi momenta 
kp as function of the laboratory energy E^. The free cross section (kp = 0) is 
compared to the experimental total np cross section [15]. 

In the present work we use both, mean fields and in-medium cross sections ob- 
tained in DB calculations of the Tubingen group [11,14] utilizing the Bonn A 
one-boson-exchange potential as the bare NN interaction V which yields very 
reasonable saturation properties for infinite nuclear matter. Fig. 1 shows the 
energy dependence of the in-medium neutron-proton (np) cross section [14] at 
Fermi momenta kp = 0.0, 1.1, 1.34, 1.7/m -1 , corresponding to p ~ 0, 0.5, 1, 2po 
(po — 0.16/m -3 is the nuclear matter saturation density). The presence of the 
medium leads to a substantial suppression of the cross section which is most 
pronounced at low laboratory energy E\ & \> and high densities where, in addi- 
tion to the (M*/M) 2 scaling of (3) the Pauli-blocking of intermediate states 
is most efficient. At larger asymptotic values of 15-20 mb are reached. 
However, not only the total cross section but also the angular distributions are 
affected by the presence of the medium. The initially highly forward-backward 
peaked np cross sections becomes much more isotropic at finite densities [14] 
which is mainly do to the Pauli suppression of soft modes (7r-exchange) and 
corresponding higher partial waves in the T-matrix. 

Heavy ion collisions are described by the relativistic Boltzmann-Uehling-Uhlenbeck 
(RBUU) transport equation, which has the form [12] 



k*"^ + {k* v F" v + M*d£M*) d k J I f(x, k*) ' 



2(2tt) 9 

x / [hhfh - fhhh] (4) 

J K 2 K 4 

where the single particle distribution function is given by fi = /(x, t; k*) and 
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the hole-distribution by fa = (1 — /(x, t; k*)). The collision integral accounts 
explicitly for final state Pauli-blocking while the in-medium scattering ampli- 
tude includes the Pauli-blocking of intermediate states. 

The dynamics of the lhs of eq.(4), the drift term, is determined by the mean 
field. Here the attractive scalar field £5 enters via the effective mass M* and 
the repulsive vector field via kinetic momenta /c* and via the field tensor 
p^v _ q^Yj v — £?"E M . The in-medium cross sections enter into the collision 
integral via the transition amplitude 

W = (2tt) 4 <5 4 (k + k 2 - k 3 - k 4 ) (M*) 4 |T| 2 (5) 
with relation (3). 

It was found that in relativistic heavy ion collisions the phase space distri- 
bution is not equilibrated during much of the reaction [4]. This modifies the 
in-medium T-matrix and thus the mean field and the cross section via the in- 
termediate state Pauli operator. To account for this approximation, the mean 
field is described in the Local (phase space) Configuration Approximation 
(LCA) where the anisotropic phase space is taken into account locally in a 2- 
Fermi-sphere configuration [4]. In a fully consistent treatment one should treat 
the in-medium NN scattering in the same way, i.e. for a scattering of a nucleon 
to a 2-Fermi sphere configuration which is not available yet. To make the dis- 
cussion more transparent here, we will discuss results of transport calculations 
where the mean field is always treated in the LCA approximation. 

The connection of the microscopic BUU approach to a hydrodynamical pic- 
ture is established via macroscopic transport coefficients such as the shear 
viscosity 77 and the heat capacity [16]. Though defined close to the equilib- 
rium limit, transport coefficients give insight in the bulk properties of excited 
nuclear matter. In particular the shear viscosity is closely connected with the 
in-medium cross section. Quantum effects enter via the final state Pauli block- 
ing for final states and for intermediate states in the BS-equation (1). The cross 
section is not only reduced at finite density but also becomes more isotropic 
which increases the average momentum transfer and reduces the viscosity of 
the system. 

It is of great interest where such effects are experimentally accessible in heavy 
ion reactions. While the close connection of collective flow observables to the 
mean field is well established [1,7,8,4], the experimental access to in-medium 
modifications of the cross section is not so straightforward. Stopping observ- 
ables are in principle natural tools to study the effectiveness of binary colli- 
sions. However, the comparison to experimental rapidity distributions up to 
now did not deliver conclusive information on in-medium modifications of the 
cross section. 
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Very recently, the FOPI collaboration at GSI performed a systematic analysis 
of nuclear stopping in the most central reactions over a wide energy range 
form 0.1 -J- 1.5 AGeV [17]. The complete phase space coverage allowed not 
only the measurement of longitudinal rapidity distributions but also those 
transverse to the beam direction. The data showed that even in most central 
reactions the system is not completely equilibrated which is reflected in a 
broader longitudinal distribution. A quantitative measure for the equilibration 
was introduced by Reisdorf et al. [17] by the ratio 

Vani = -aHz) (6) 



of the longitudinal (z) and transverse (y) variances of the rapidity distribu- 
tions. Hence, a value of var t i < 1 indicates incomplete stopping. It was further 
observed that stopping (in central reactions) and transverse flow (in semi- 
peripheral reactions) show a strong correlation in that the excitation function 
for both observables has a broad maximum around 0.4-0.6 AGeV [17]. 

Fig. 2 shows the incident energy dependence of the observable var t \ and the 
scaled directed flow P^°J ir = sign(yi)p^) /P pro j, where the sum runs over all 
charged particles. A comparison to experimental data [17] requires an accurate 
determination of the centrality class in the same way as in the experiment. In 
particular, this is important for the observable var t \ due to its high sensitiv- 
ity to the centrality selection. We have performed a systematic study on the 
impact parameter dependence of the flow and stopping observables. -P^L was 
theoretically determined both for the centrality interval PM4 (top-right panel 
in Fig. 2), corresponding to b=2-6 fm and for the fixed impact parameter b=4 
fm (top-left panel in Fig. 2), corresponding to the maximum magnitude of di- 
rected flow. The last condition was used in determining P^ ir experimentally 
[17]. A consistent determination of the stopping observable v ar\\ and the high- 
est centrality class is more difficult due to its high sensitivity on the centrality 
selection. In Ref. [17] the observable ERAT was used for the experimental 
determination of central collisions. Comparisons with IQMD simulations [18] 
provide an impact parameter range b ~ 2 — 2.5 fm, consistent with previous 
studies [4]. However, because of numerical fluctuations of the ERAT distribu- 
tion a determination of the impact parameter to better than 5b = 0.25 — 0.5 
fm is not possible, which essentially affects the final result of var t \. The fluctu- 
ations of var t \ according to these Sb are indicated in Fig. 2 (bottom) as error 
bars. Furthermore, as for the experimental data [17], var t \ is determined for 
central reaction (b < 2.5 fm) by Gaussian fits to the final rapidity distribu- 
tions. The theoretical calculations are performed for two cases, i.e. using free 
and in-medium total and differential NN cross sections (solid and dashed lines 
in Fig. 2, respectively). 

Generally the microscopic calculations reproduce the qualitative behavior of 
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Fig. 2. Energy dependence of the scaled directed flow P^J ir = Pxdir/Ppxo] (top pan- 
els) and the degree of stopping var t \ (bottom) for Au+Au collisions. The directed 
flow and the observable varti are experimentally determined for the scaled impact 
parameter = b/b ma x ~ 0.4 (b max being the maximal value of the impact pa- 

1/3 1/3 

rameter defined as bfyiax — l*l^(^proj ^targ)) corresponding to the magnitude of 
maximal global flow and for b^ ' < 0.15 in semi-central and central Au+Au colli- 
sions, respectively. Transport calculations using free (solid lines) and in-medium NN 
cross sections (dashed lines) are compared to data from the FOPI collaboration [17] 
(filled diamonds). The directed flow is theoretically determined for a fixed impact 
parameter of b = 4 fm (top-left panel) and for the centrality class PM4 [b = [2, 6] fm) 
(top-right panel). The observable varti (bottom panel) is theoretically calculated for 
b < 2.5 fm with a variation of 5b = 0.5 fm in the impact parameter range. 



the data. The stopping power as well as the scaled directed flow have a broad 
maximum as a function of beam energy. The use of in-medium cross sections 
has only little impact on the collective flow as expected. The flow is slightly 
reduced but at high energies the data are still somewhat overestimated. This 
is understood from the momentum dependence of the nucleon-nucleus optical 
potential. Within the microscopic Dirac-Brueckner approach [11] the optical 
potential is too repulsive at energies above 0.6 AGeV and does not reproduce 
the empirically observed saturation of the optical potential [4]. In heavy ion 
collisions this leads to an overestimation of the collective flow, which has 
been cured adopting ad-hoc phenomenological parameterizations of the optical 
potential [8]. The main reason for the decrease of the flow at high energies 
is the opening of inelastic channels (pion, kaon production etc.). Such in- 
elasticities are included in the collision part of the transport calculation but 
not on the level of the mean field. The theoretical understanding of the bare 
NN interaction above the inelastic threshold is still an open issue [19]. 
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Fig. 3. Elastic and inelastic NN cross section multiplied by average Pauli blocking 
factors for the initial phase space configurations of heavy ion reactions. 

The energy dependence of varw can be understood more easily in a qualitative 
way: at low energies Pauli blocking is most efficient, suppressing 2-body colli- 
sions which leads to a higher transparency. With increasing energy the block- 
ing becomes ineffective whereas the elementary cross section drops. Hence the 
maximum reflects the combination of final state Pauli blocking and a decreas- 
ing cross section. For a quantitative description of the data the additional 
reduction of the in-medium cross section by Pauli blocking of intermediate 
states and M* scaling is seen to be essential. 

The interplay between Pauli blocking and the NN cross section can be seen 
from Fig. 3. This figure shows the isospin averaged free cross section multiplied 
by average Pauli blocking factors for the initial local momentum configurations 
in the central cell of a heavy ion reaction. These are given by two nuclear mat- 
ter currents of density po separated by the beam velocity which in momentum 
space correspond to Lorentz boosted Fermi ellipsoids. The blocking factors are 
evaluated for two representative nucleons in the centers of the ellipsoids. One 
clearly sees how the interplay between Pauli blocking and decreasing cross 
section leads to a broad maximum in the effective cross sections at about 0.5 
GeV which appears in the same region as that of vara- 

It is also important to obtain an understanding of the stopping behavior in 
terms of nuclear matter bulk properties. However, ideal one-fluid hydrodynam- 
ics which is most directly related to nuclear matter bulk properties completely 
fails at SIS energies [17]. In contrast, quantum effects such as Pauli blocking 
of two-body collisions are essential for the degree of stopping. In Ref. [16] 
nuclear transport coefficients have been derived from the microscopic BUU 
equation. In the Boltzmann statistical limit these correspond to the first or- 
der Chapman-Enskog coefficients. At large density the microscopic transport 
coefficients differ, however, essentially from the Chapman-Enskog results [16]. 
According to [16] the shear viscosity of a nuclear liquid can be parameterized 
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Fig. 4. Shear viscosity of nuclear matter as a function of density p/po (indicated 
by labels) and temperature. Diamonds show the values obtained in central Au+Au 
reactions at the stage of maximal compression. The values correspond to incident 
energies of 0.15, 0.25, 0.4, 0.6 and 0.8 AGeV (with rising T). 

as a function of p and T 

1700 (p\ 2 22 (P\°\ 5.8y/T 

V ~ T 2 [po) l + T-10- 3 ^po/ 1 + 160T- 2 U 



Fig. 4 shows 7] according to eq. (7) as a function of T for various p together 
with the values which have been extracted from central Au+Au reactions. The 
T — p for which rj is evaluated have been determined in the central cell at the 
stage of maximal compression. The temperature was thus extracted by fits of 
Fermi distribution to the local phase space as delivered from the transport 
calculations [3]. From Fig. 4 one sees that rj exhibits a minimum as function 
of beam energy at E'Lab = 0.25 AGeV which is close to the maximum of var t \. 

In summary, we performed for the first time RBUU calculations where the 
mean field and the in-medium cross section are consistently based on the same 
in-medium Dirac-Brueckner T-matrix. We investigated the recently proposed 
observable var t \ as a measure for the degree of stopping and equilibration in 
central reactions. In-medium effects on the NN cross section are essential in 
order to reproduce these data. Qualitatively, the behavior of the var t \ excita- 
tion function can be understood by the interplay between the decreasing NN 
cross section and final state Pauli blocking. The interpretation of nuclear bulk 
properties in terms of transport coefficients indicates that the shear viscosity 
is minimal in the regime where maximal stopping occurs. 
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